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Abstract: VERITAS employs a 12m segmented mirror and pixellated photomultiplier tube camera to 
detect the brief pulse of Cherenkov radiation produced by the extensive air shower initiated by a cosmic 
high-energy gamma ray. The VERITAS data acquisition system consists of a 500 Mega-Sample-Per- 
Second custom-built flash ADC system, which samples the Cherenkov light pulse every 2 nanoseconds. 
The integrated charge in each flash ADC channel is proportional to the amount of Cherenkov light incident 
on the corresponding photomultiplier tube. Accurate reconstruction of the integrated charge is required 
for accurate energy estimation and spectral reconstruction. A reliable calculation of the integrated charge 
at low intensities can lead to a reduction in the energy threshold of the system, and an increase in sensi- 
tivity. This paper investigates and compares several approaches for evaluating the integrated charge. The 
Cherenkov pulse timing information in the flash ADC readout has the potential to assist in background 
rejection techniques. Various methods for extracting the timing information are investigated and excellent 
timing resolution is achieved. 



' Introduction 

, There are many methods [6] which can be used to 
• evaluate the digitised Cherenkov signal produced 
■ by the VERITAS FADC [2] system. In this paper, 
five such methods (referred to as trace evaluators) 
are described and the charge integration character- 
istics of each method compared using a Monte- 
Carlo simulated photon data set. This study could 
' aid accurate reconstruction of low-intensity events, 
' which is one of the most challenging aspects of the 
analysis of Cherenkov telescope data. Laser [5] 
calibration data are used to compare the timing res- 
olution inherent to each trace evaluator, and a dig- 
ital processing scheme which can further enhance 
the timing resolution is introduced. These methods 
have been developed and implemented with VE- 
GAS [3]. 

Methods 

In this section each trace evaluator will be de- 
scribed and the manner in which the integrated 
charge and pulse arrival time is calculated is dis- 



cussed. The integrated charge is defined as the sum 
of the trace in digital counts over some integration 
window. The pulse arrival time (hereafter Tq) is 
defined as the time at which the pulse reaches 50% 
of its absolute maximum. 

The first method is the simple-window trace evalu- 
ator which assumes a-priori knowledge of the loca- 
tion of the Cherenkov pulse in the readout window. 
The second method is the dynamic-window trace 
evaluator which improves on charge integration by 
sliding an integration window along the readout 
window to seek the Cherenkov pulse. The first two 
evaluators only calculate Tq to the nearest sam- 
ple. The third method is the linear-interpolation 
trace evaluator This is not significantly different 
in terms of charge integration, but substantially im- 
proves on the calculation of Tq. The fourth method 
is the trace-fit trace evaluator which fits the follow- 
ing function to each trace: 

{qo exp for t<ta 

go cxp for t>to 

(1) 
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In this equation, q{t) is the FADC charge at time t, 
to is the time of trace maximum, qo is the trace 
amplitude at t — to, and a and a are param- 
eters describing the shape of the trace. This fit 
function essentially has an asymmetric-Gaussian 
shape and improves the calculation of To over the 
simple-window method. The fifth method is the 
matched-filter trace evaluator which uses a digi- 
tal filter based on the assumed shape of the FADC 
pulse to integrate the charge. The matched-filter 
trace evaluator is a somewhat more sophisticated 
than the other methods, thus it is described here in 
more detail. 

A matched filter is so called because its shape is 
defined by the expected form of the received data. 
The matched filter's pulse shape is a time-reversed 
version of the expected pulse shape. Thus for an 
expected pulse shape h{t), the ideal matched-filter 
hm{t) is 

K,{t) = h{T-t) (2) 

for < i < T where T corresponds to the end of 
the trace. The output from a filtering application 
is calculated by a convolution of the input with the 
filter 

y{t)^ [ r{t)h„,{T - t) dt (3) 
Jo 

which reduces to the cross correlation of r{t) and 
h{t) with zero lag. 

y{t)= ( r{t)h{t)dt (4) 

In order to construct the matched filter a stan- 
dard laser calibration run is used which is normally 
used to flat field the camera. For each event, and 
for each channel, a section of the laser pulse is 
extracted, and aligned relative to some predeter- 
mined point. This extracted pulse is summed for all 
events for each channel. The summed trace is nor- 
maUsed (Figure 1), and Fourier transformed. The 
filter is applied to the FADC data by multiplying 
the Fourier transform of the FADC trace (denoted 
S{lli)) with the conjugate of the filter transform, 
H{lo), and then applying an inverse Fourier trans- 
form 
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Figure 1 : Normalised subset of FADC trace used 
to construct a matched filter. 



y{t) = T-^ [S(uj) X H{uj)\ (5) 

which yields the cross correlation function y{t). 
The maximum of the cross correlation is propor- 
tional to the integrated charge of the FADC trace. 
In order to establish the constant, a series of special 
laser calibration runs is taken with continuously in- 
creasing laser intensity. The integrated charge as 
measured using the dynamic -window trace evalua- 
tor is compared to the output of the matched-filter 
trace evaluator and used to establish the constant 
[4]. 

When analysing data, the charge from a trace is 
derived by applying the matched-filter trace evalu- 
ator, and multiplying the output by the appropriate 
constant for that channel. The pulse arrival time 
is determined by the location of the maximum of 
the cross correlation, thus the arrival time can only 
be determined to the nearest FADC sample (much 
like the simple-window ti-ace evaluator). 

Integral Charge Evaluation 

In order to examine the charge evaluation quality 
of each trace evaluator, a data set of photon impacts 
on the camera is simulated. The arrival time of the 
photons is assumed to be Gaussian with an RMS of 
2 ns. The simulation is performed using GrISUDet 
[1]. A comparison of an FADC trace simulated in 
this way with a real trace from a laser calibration 
run is shown in Figure 2. 
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Figure 2: Comparison of real and simulated FADC 
trace. The real trace is indicated by the points. 

The simulated data set is divided into subsets such 
that each subset only has events with a certain 
number of photoelectrons. Data sets with from 
one to thirty photoelectrons are generated in this 
way. This allows the charge reconstruction as a 
function of the known number of photoelectrons 
to be evaluated. For each trace evaluator, a distri- 
bution of integrated charges (in digital counts) is 
generated for each photoelectron multiplicity. In 
terms of charge evaluation, the quality of the trace 
evaluator is determined by the RMS of the distri- 
bution of integrated charges for a constant input. 
The difference in the RMS of the simple-window 
trace evaluator, and each other trace evaluator as a 
function of the number of photoelectrons is shown 
in Figure 3. Thus, the RMS of the simple-window 
trace evaluator is used as a baseline against which 
the other trace evaluators can be compared. For 
small pulses ( < 5 photoelectrons), the matched- 
filter trace evaluator provides the smallest RMS, 
however the RMS quickly increases with the num- 
ber of photoelectrons. This is to be expected as 
although small pulses are dominated by noise, the 
matched filter is is able to pick out the signal from 
the trace. Conversely, the trace-fit trace evalua- 
tor gives a very large RMS for small pulses. This 
is attributed to ill-fitting of small, poorly-defined 
pulses. At approximately four photoelectrons, all 
the trace evaluators yield similar results. Beyond 
that, the trace-fit trace evaluator is superior, and 
only the matched-filter trace evaluator is signifi- 
cantly worse. 



Figure 3: Comparison of charge resolution relative 
to the simple-window trace evaluator 
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Figure 4: Comparison of original and resampled 
FADC trace. 

Trace Resampling 

One tool commonly used in digital signal process- 
ing is resampling in the time domain. The resam- 
pling is achieved by applying a Fourier transform 
to the FADC trace, zero-padding in the frequency 
domain, and applying an inverse Fourier transform. 
Zero padding in the frequency domain is achieved 
by simply adding zeros to the end of the Fourier- 
transformed trace. This has the effect of setting 
higher frequencies to have zero amplitude. The 
inverse Fourier transform results in a trace which 
has been resampled in the time domain. Figure 4 
shows a comparison between a raw and resampled 
trace. 
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Figure 5: Timing resolution as a function of FADC 
trace size. 

Timing Resolution 

The timing resolution is determined by how well 
the arrival time of an asynchronous laser flash in- 
cident on the camera plane can be measured. 

The timing resolution is defined as the width of 
a Gaussian function fit to the distribution of mea- 
sured differences between event arrival and chan- 
nel arrival time for a series of laser pulses for each 
channel. The event arrival time is defined as the 
average arrival time of all the channels in the cam- 
era. The timing resolution as a function of inte- 
grated charge is shown for all methods in Figure 5. 
As expected, the timing resolution improves as a 
function of pulse size, as for larger pulses the time 
structure is better defined. Above an integrated 
charge of 30 dc, the timing resolution is roughly 
linear as a function of charge. The timing resolu- 
tion for pulses greater than 30 dc is shown in Table 
1. 



As expected, the linear-interpolation trace evalu- 
ator has excellent timing resolution for all trace 
sizes. The matched-filter trace evaluator is excel- 
lent for small traces, however as pulse arrival times 
can only be calculated to the nearest FADC sample, 
it is not as good for large pulses. The trace-fit trace 
evaluator has poor timing resolution for very small 
pulses - mirroring the effect seen with the study 
of charge resolution, indicating that the fit func- 
tion is not suited to small pulses. For large pulses 
the trace-fit evaluator has a superior resolution. 
However, the best timing resolution is achieved us- 
ing a combination of the resampling technique and 
the linear-interpolation trace evaluator Together, 
these fast robust techniques provide a timing reso- 
lution of just 0.22 ns with these data. 

Conclusions 

Five trace evaluation techniques and a digital fil- 
tering technique have been described and com- 
pared using real and simulated data. The matched- 
filter technique holds promise for the evaluation of 
small sub-threshold events. Optimal timing resolu- 
tion has been achieved using an FADC resampling 
technique in concert with the linear-interpolation 
trace evaluator. 
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Evaluator 


Resolution (ns) 


Error (ps) 


Simple Window 


0.77 


3.2 


Linear Interpolation 


0.45 


2.9 


Trace Fit 


0.46 


2.9 


Matched Filter 


0.91 


2.7 


Resampling 


0.22 


10.9 



Table 1: 
30 dc. 



Timing resolution for pulses greater than 
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